Submitted April 7, 2000; Accepted June 28, 2000; to appear in ApJ Dec. 1, 2000 
Preprint typeset using lAT^^X style cmulatcapj v. 04/03/99 



the X-RAY SIZE-TEMPERATURE RELATION FOR INTERMEDIATE REDSHIFT GALAXY 

CLUSTERS 

Joseph J. Mohr^-^, Erik D. Reese^, E. Ellingson^, Aaron D. Lewis^ & 

August E. Evrard"* 
Submitted April 7, 2000; Accepted June 28, 2000; to appear in ApJ Dec. 1, 2000 



o 
o 
o 

(N 



(N 
> 
(N 

O 
O 

o 



Oh- 
I ' 

O 



X 



ABSTRACT 

We present the first measurements of the X-ray size-temperature (ST) relation in intermediate redshift 
{z ~ 0.30) galaxy clusters. We interpret the local ST relation (z ^ 0.06) in terms of underlying scaling 
relations in the cluster dark matter properties, and then we use standard models for the redshift evolution 
of those dark matter properties to show that the ST relation does not evolve with redshift. We then use 
ROSAT HRI observations of 11 clusters to examine the intermediate redshift ST relation; for currently 
favored cosmological parameters, the intermediate redshift ST relation is consistent with that of local 
clusters. Finally, we use the ST relation and our evolution model to measure angular diameter distances; 
with these 11 distances we evaluate constraints on VLm and f2A which are consistent with those derived 
from studies of Type la supernovae. The data rule out a model with = 1 and JIa ~ with 2.5f7 
confidence. When limited to models where JIm + f^A = 1, these data are inconsistent with VIm — 1 with 
3cr confidence. 

Subject headings: galaxies: clusters: general — intergalactic medium — cosmology 



1. INTRODUCTION 

Nearby galaxy clusters exhibit a tight correlation be- 
tween X-ray isophotal size and em ission weighted intra- 
cluster medium (ICM) temperature ( Mohr fc Evrard 1997 ; 
hereafter ME97). This correlation is evidence of regular- 
i ty; it exists in an X-ray flux limited sample of 45 clusters 
( Edge et al. 1990| ), where no attempt has been made to use 
the X-ray morphologies to exclude clusters showing signa- 
tures of recent, major mergers. The scatter around the X- 



'ay size temperature (ST) relation is approximately 159l 
in oizo, comparable to the ocattor of elliptical and lonti e- 
ular galaxies around thei r fundamental plane (J0rgensen, 
Franx fc Kjargaard 1996| ) . This small scatter in the galaxy 



cluster scaling relation is intriguing, because (1) there is 
overwhelming evidence that galaxy clusters are still ac cret- 
ing mass (e.g. [Mohr et al. 1995| , |Buote fc Tsai 1996| ) and 
(2) elliptical galaxies are generally thought to be among 
the most regular objects in the universe. 

ME97 use 48 N-body and hydrodynamical simulations 
of cluster formation in four different cosmological models 
to address this apparent contradiction between regularity 
and ongoing accretion in nearby clusters. Using simula- 
tions from both Q.m = 0.3 and VIm = 1 cosmologies, they 
show that a tight ST relation is expected even in cosmogo- 
nies where there is significant cluster growth at the present 



ipoch. 



The high dcgrGC of regularity implied by the ST relation 
io Durprioing, bocauGO the well known correlation between 



X-ray luminosity and emission weighted mean te mpera- 
ture (the Lx — Tx relation) has very large scatter ( Davicj 
et al. 1993 ). ME97 show that the same cluster ensemble 
15% scatter in the ST relation exhibits 
Tx relation. This 



which exhibits a 
a 52% scatter in Lx around the Lx 



higher scatter in the Lx — Tx relation results from the 
sensitivity of the X-ray luminosity to the densest regions 
of the cluster- a sensit ivity to the presen ce or absence of 
so-called cooling flows ( Fabian et al. 1994 ) . This interpre- 
tation is supported by more recent work where cluster en- 
sembles specially chosen to contain no cooling flow clusters 
conform to Lx — Tx relations with signi ficantly reduced 
scatter of 25% ( [Arnaud fc Evrard 1999| ). Additionally, 
when central parts of cooling flow clusters a re excluded, 
the s catter in the Lx — Tx relation decreases ( Markevitch 



1998) 



Finally, cluster regularity is also evident in the tight cor- 
relation between ICM mass and emission weighted temper- 
ature (the MjcM — Tx relation). When measuring ICM 
mass within a limiting radius of rsoo (the radius where 
the enclosed overdensity is 500 times the critical density) 
the scatter in mass about the Mjcm — Tx relation is 17% 



( Mohr, Mathiesen fc Evrard 1999 ). Observational studies 
of the ST relation followed by work on the Lx — Tx and 
Mic M — Tx relations support a scenario where clusters ex- 
hibit regularity similar to that of elliptical galaxies in the 
properties measured on the scales of their virial regions, 
but exhibit significant irregularities in the properties of 
the densest, central regions where physical processes other 
than gravity and gas dynamics- such as radiat ive cool- 
ing and magnetic fields- play significant roles (Mohr fc 



Evrard 1991, Arnaud fc Evrard 1999, Mohr, Mathiesen 



Evrard 1999[). The evidence indicating cluster regularity 



is balanced by evidence for departures from regularity; the 
scatter in the observed scaling relations is larger than can 
be accounted for by the measurement uncertainties. This 
resolved scatter contains clues about, among other things, 
cluster peculiar velocities and departures from equilib- 



rium. 
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X-RAY ST RELATION AT INTERMEDIATE REDSHIFT 



Here we examine the ST relation at intermediate red- 
shift (0.19 < z < 0.55) using ROSAT HRI observations 
of the Canad ian Network for Observational Co smology 
(CNOC; e.g. |Yee et al. l"996| [Lewis et al. 1999D cluster 
sample. The ST relation provides a potentially powerful 
tool to study the expansion history of the universe. As ex- 
plained in detail below, the ST relation is rather insensitive 
to cosmological evolution. Thus, armed with an accurate 
model of cluster evolution, one could use the ST relation 
to me asure distances at intermediate redshift, constraining 



the dc celeration parameter qq. 



We first present the X-ray ST relation in the nearby 
cluster sample (§|^) and then present an interpretation 
of the ST relation in terms of regularity in the under- 
lying dark matter properties of the cluster. Section || 
describes observatio ns of the ST relation in intermediate 
redshift clusters. In S3.3 we use these observations to con- 



strain cosmological parameters. Section ^ contains a sum- 
mary of our conclusions. Throughout the paper we use 
Hq = 50/i5okm/s/Mpc. 
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Fig. 1. — The X-ray isophotal size versus emission weighted mean 
ICM temperature Tx for an X-ray flux limited sample of 45 nearby 
clusters. We use the isophote I = 3.0 X 10"^"* erg/s/cm^/arcmin^ 
within the cluster rest frame 0.5:2.0 keV band. The solid line rep- 
resents the best fit relation, and the RMS scatter about this line is 
15% in size. 

2. LOCAL X-RAY ST RELATION 

Below we present observations of the ST relation in an 
X-ray flux limited sample of clusters and then discuss clus- 
ter scaling relations. 

2.1. Observations 

Fig |l contains a plot of the X-ray ST relation for an 
ensemble of nearby galaxy clusters (see ME97 for clus- 
ter list). This sample is essentially a low redshift sample, 
Zmedian ^ (z) = 0.055 with RMS variation cr^ = 0.03, but 
the full redshift range is 0.01 < z < 0.19. The cluster 
isophotal size Rj is plotted versus the emission weighted 
mean ICM temperature Tx, where 



and Aj is the area enclosed by the isophote /. The 
isophotal size is a model-independent measure extracted 
directly from the X-ray image. ICM temperatures come 
from the literature (listed in ME97). In this figure we 
use the galactic absorption corrected surface brightness 
/ — 3.0 X lO^^^'ergs s~^cm~^arcmin~^ within the clus- 
ter rest frame 0.5-2.0 keV band. In the conversion from 
count rate to physical flux units we use PROS, assume a 
Raymond-Smith spectrum with the publis hed mean tem- 



perature Tx and ^ cosmic abundances (Mushotzky & 



Loewenstein 1997). We convert between the angular size 
and physical size of the isophotal region using the angular 
diameter distance, and we also correct / for cosmological 
dimming (1 -|- z)^. 

The best fit power law to this local ST relation has the 
form 



Ri = (0.71 ±0.02) 



T^\ 
6keVy 



1.02±0.11 



hj Mpc 



(2) 



and the scatter about this relation is 15% in Rj. This fit 
is shown as a solid line in Fig |^. Because this is the local 
ST relation, it is only weakly dependent on the deceler- 
ation parameter. For this relation we have assumed the 
parameters ftj^j = 0.3 and JIa = 0.7; for il^ = 0.3 and 
fix = the normalization is 1.8% lower, and for CIm ~ 1 
the normalization is 2.7% lower. 

2.2. Cluster Scaling Relations 

Galaxy cluster scaling relations between, for example, 
virial mass and temperature T are expected if clusters 
are approximately virialized. The spherical collapse model 
predicts that newly collapsed objects will have mean den- 
sities which are A times the critical density pcrit = 
3H^/8ttG, where H is the Hubble parameter and in an 
Einstein-de Sitter model A = IStt^. Therefore, we define 
the virial radius ta of an object with virial mass Ma to 
be 



''A = 



Ma 



^nApcrit 



1/3 



(3) 



Newly collapsed clusters may also satisfy the virial relation 
GMa = aTrA (4) 

where T is the virial temperature and a is a number de- 
pendent on the cluster density and temperature structure; 
in the case that cluster structure is self similar (i.e. that 
massive clusters are simply rescaled versions of low mass 
clusters) a will be constant for all T and Eqns |^ and ^ lead 
to the well known scaling relations between virial mass, ra- 
dius and temperature 



Ma = 



and TA = 



y/Apcr 



(5) 



Ri - ^Ai/tt 



(1) 



where a' is a dimensionless structure parameter like a. 
These scaling relations are expected to hold for a range 
of rs and Ms, where rs is the region with enclosed over- 
density S with respect to pcrit and Ms is the corresponding 
enclosed mass. Numerical cluster simulations indicate that 
the emission weighted mean ICM temperature Tx is an ad- 
equate proxy for the virial temperature T, and that these 
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scaling relations exist even in cluste r populations which 
exhib it merging and substructure ( Evrard, Metzler & 
Navai ro 1996 , Bryan fc Norman 1998 ). Using simulations 
to normalize these relations, one finds that clusters with 
Tx of 10 keV have virial masses M200 = 4 x lO^^h^^MQ 
and virial radii r2oo = 4:h'^Q Mpc. The small scatter in the 
observed ST relation, the correlation between ICM mass 
and temperature (sec Section 1), and a direct study of the 
correlation between binding mass (estimated using clus- 
ter galaxy kinematics) and emission weighted temperature 
( Horner, Mushotzky fc Scharf 1999 ) provide evidence that 
the conditions required for Eqn |5| to be valid are generally 
met in real galaxy clusters. 

2.3. X-ray ST Relation 

The X-ray surface brightness (units: ergs s~^cm~^sr^^) 
at a projected radius Rs from the cluster center is 



HRs) 



27r(l + z)-" Jo 



dlne{r)nH{r)K{T) (6) 



where rif. and jifj are the electron and proton number den- 
sities, r = R^g is the distance from the cluster center 
and A(r) is the X-ray emission coefficient which includes 
contr ibutions from thermal bre msstrahlung and line emis- 
sion (Raymond & Smith 1977). We can express the ICM 



densit y pg in terms of the underlying dark matter density 
Pdm'- Pg[r) = tgPdm{r)g[r), where /g is the 1(JM mass 
fraction within the virial region and g{r) is a function de- 
scribing differences in the dark matter and ICM density 



distril )utions; then, using — Pg/ pifUp, we can rewrite 
the X-ray surface brightness as 



I{Rs 



2'Km'^p.ePH (1 



dlpl„Xr)9\r) (7) 



We bring the emission coefficient out of the integral and 
use tile value at T — Tx', although this is only strictly 



valid for an isothermal gas, the temperature msensitivity 
of A(r) band limited to 0.5:2.0 keV makes this an ex- 
cellent approximation eve n in the presence of departures 
from isothermality (e.g. Fabricant, Lecar fc Gorensteirj 



1980i |Mohr, Mathicscn fc Evrard 1999| y 

Finally, we express the dark matter density profile in 
terms of the characteristic overdensity of the virial re- 
gion: pdm{r) = Apc„tfdm{rAy), where fdmiv) encodes 
the dependence of the dark matter profile on a dimen- 
sionless radius y = r/r^; this approach is consistent with 
numerical cluster simulations, which indicate that clus- 
ter dark matter density profiles have a "universal" form 
dNavarro, Frenk fc White 1997| , hereafter NFW). More 
recently, higher resolution simulations have shown differ- 
ences between the inner profile (at radii which are 1% of 
the virial radius) and the form proposed by NFW ( Moore 
et al. 1999 , Jing fc Suto 1999| ); these differences appear to 
be smallest for cluster scale haloes, and none of our con- 
clusions are sensitive to the behavior of the density profile 
at 1% of the virial radius. Expressing p^m in this way, and 
removing the characteristic scale Rs from the integral, the 
X-ray surface brightness becomes 



HRs) = 



f^AiT)pl„,A^Rs 

2'Kmlp,eP.H (1 + zf 



e 



(8) 



where is a dimensionless integral which encodes the 
shape of the ICM density profile: 6 = d\ f2„-,{ri)g'^{rj), 

where 77 = Vl + and = {r/Rsf - 1. 

We relate Rs and Rj using the shape of the typical X- 
ray surface brightness profile; I{R) is typically well fit b y 
the so-called (3 model (|Cavaliere fc Fusco-Femiano 197S) 



I{R) = lo 




-3/3+1/2 



(9) 



where Rc is the core radius; well outside the core the sur- 
face brightness falls as I{R) cx R^~^^. Therefore we write 



Ri = Rs 



f HRs) \ 



1/(6/3-1) 



(10) 



Using Eqn. ^ to substitute Tx for Rs, one can readily de- 
termine the ST relation at a particular redshift in the case 
that cluster structure is self similar 



Ri cx T^ where a 



3f3 



6/3-1 



(11) 



For /3 = 2/3, a typical observed value (Jones fc Forman 



1984| , |Mohr, Mathiesen fc Evrard 19991 ), the predicted ST 
relation has a = 2/3, shallower than the observed rela- 
tion: a = 1.02 ±0.11. The observed increase in ICM mass 



fraction /q with cluster temperature Tx (Mohr, Mathiesen 



fc Evrard 1999) is enough to explain the steepness of the 
observed relation. The physics underlying this systematic 
variation in cluster gas mass fraction is likely heating of 
the intergalactic medium during the process of galaxy for- 
mation. This so-called preheating affects the structure of 
the gas in low mass clusters more so than in high mass 
clusters Je_ 



fc Tozzi 199S , 



Metzler fc Evrard 1994, Cavaliere, Menci 



onman. Cannon fc Navarro 199S| ). As dis 



cussed below, in comparing the local and intermediate red- 
shift ST relations it is important only that the effects of 
this preheating be similar in both samples. 

3. INTERMEDIATE REDSHIFT X-RAY ST RELATION 

Before presenting the intermediate redshift ST relation, 
we present an evolution model and our measurement meth- 
ods. 



3.1. Evolution of the ST Relation 

Comparison of ST relations in nearby and distant pop- 
ulations requires a model for how the expected changes in 
cluster structure with increasing redshift will affect the ST 
relation. In the previous section we presented an explana- 
tion of the correlation between cluster X-ray isophotal size 
and emission weighted ICM temperature Tx as a manifes- 
tation of underlying scaling relations in the dark matter 
properties of the cluster. Following this line of reasoning, 
we now express the evolution of the ST relation in terms 
of the redshift evolution of these same dark matter prop- 
erties. 

The evolution of the dark matter scaling relations is 
apparent in Eqn ^ The normalization of these scaling re- 
lations changes with redshift according to the change in 
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Pcrit and the change in the characteristic overdensity A of 
collapsed haloes. Within Einstein-de Sitter models A does 
not evolve, whereas evolution is expected in low JIm mod- 
els. For simplicity of presentation, in the following analysis 
we explicitly follow only changes in pcriti but accounting 
for changes in A would not affect our conclusions. Using 
H{z) = HqE{z), we express the evolution of the mass and 
size scaling relations as 

Ma(T,0) , , RAiT,0) , , 

MA(r,z)^ ^ andi?A(T,z) = t!, I (12) 



E{z) 



E{z) 



where E'^{z) = VLm^^ + zf -h (1 -VLm - f^A)(l + zf -I- J^a, 
and and J^a are the present epoch contributions to 
the density parameter from matter and "dark energy" , re- 
spectively. These simple evolution models are valid only 
for self similar evolution- i.e. in the case that distant clus- 
ter dark matter profiles are structurally similar to those of 
nearby clusters. 
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Fig. 2. — The 1, 2 and 3 u joint confidence regions for the /3 
model fit to MS-1455 (Ax^ = 2.3, 6.2, 11.8). Overplotted are fines 
of constant isophotal size Qj in arcseconds. Note that these fines 
run at shafiow angles with respect to the /3-9c correlation, reducing 
the correlation's effect on our 9/ uncertainties. 

From Eqn |^ it follows that the X-ray surface brightness 
at i?5 evolves as I{Rs) (x E^{z) and that the normalization 
of the ST relation evolves as 



Ri{T,z) = Ri{T,0) E'^{z)wheTer] 



4-6/3 
6/3-1 



(13) 



Interes tingly, for the most common cluster profile — 2/3 
(local: ^ ones fc Forman 1984| , |Mohr, Mathiesen fc Evrard| 
1999[|this intermediate redshift sample: {/3) = 0.63±0.04), 
so ry = 0, corresponding to no evolution in the ST rela- 
tion. In essence, a cluster of a given Tx is denser in the 
past, boosting its X-ray surface brightness and tending to 
increase the isophotal radius Ri, but this cluster is also 
smaller in the past, tending to decrease the isophotal ra- 
dius Rj. If the ICM density profile behaves as pg oc r~^, 
these two effects cancel. Thus, a population of clusters 



with measured /3 ~ 2/3 will exhibit an ST relation which 
is relatively insensitive to cosmological evolution. In prin- 
ciple, this behavior makes the ST relation ideal for mea- 
suring angular diameter distances to high redshift clusters. 
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Fig. 3.— The X-ray ST relation for 11 members of the CNOC 
cluster sample with measured Tx- For this figure the conversion 
from measured 0j Ri is done assuming Qm = 0.3 and Q\ = 0.7. 
Error bars denote 1 cr uncertainties in emission weighted mean tem- 
perature Tx and in isophotal size Rj. The solid line is the best 
fit ST relation determined from the nearby cluster sample (Fig hi), 
and the dashed line is the best fit ST relation for the intermediate 
redshift sample. 

3.2. Observations 

We use archival ROSAT HRI observations of 11 of the 
14 CNOC clusters discussed in Lewis et al. (1999); these 
11 clusters are those with published Tx measurements. A 
detailed description of the reduction to 0.5" pixel images 
is given in Lewis et al. (1999). The HRI angular res- 
olution is ~ 5" FWHM; we further bin these images to 
a final pixel size of 2" without significant loss of angular 
resolution. Despite the far higher median exposure time 
(33 ks compared to 8 ks), the image quality of these inter- 
mediate redshift clusters is far lower than for our nearby 
sample (due to combination of lower observer frame sur- 
face brightness and higher instrument background). The 
image quality is too poor to allow isophotal sizes mea- 
surements directly from the images using the approach 
described in S0; therefore, we fit circular (3 models (same 
form as Eqn pTwith 6 substituted for R) to these images 
using the software develope d for measuring SZE+X-ray 
distances (Reese et al. 2000). Essentially, we find the set 
of parameters /q, 0c, /3, cluster centroid {a, 6) and local 
background Ihkg which maximizes the likelihood of consis- 
tency between model and data. In all these fits we fix the 
background to the value measured in an annulus extending 
from 4.5' to 5', a background dominated region with an es- 
sentially negligible contribution from the cluster. Fits are 
performed to the central portion {0 < 4.2') of the image 
which includes the cluster and a local background region. 
In Table || we list the best fit parameters with estimates 
of the statistical uncertainties. The central surface bright- 
ness and measured background are both given in detector 
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units of cts s~^arcmin^^. 

These cluster parameters are in reasonably good agree- 
ment with those presented in Lewis et al. (1999), with 
the exception of two clusters A2390 and MS1358. The 
differ ent best fit parameters in these two clusters stem 



from ^he treatment of emission excesses; we fit the surface 
brightness model to the entire dataset, whereas Lewis et 
al. (1999) exclude the central region in both these clusters. 
This exclusion approach typically leads to larger core radii 
and correspondingly higher /3's than would fitting to the 
entire cluster as in our method. 
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Fig. 4. — The X-ray ST relation for 11 members of the CNOC clus- 
ter sample. For this figure the conversion from measured 9j Ri 
is done assuming Hjv/ = 1. Error bars denote 1 a uncertainties in 
emission weighted mean temperature Tx and in isophotal size Rj. 
The solid line is theJpest fit ST relation determined from the nearby 
cluster sample (Fig hi), and the dashed line is the best fit ST relation 
for the 11 intermediate redshift clusters. 

An examination of the fit residuals in our sample pro- 
vides no clear evidence for systematic differences between 
the data and the models with the exception of A2390, 
which exhibits an asymmetry o r centroid variation (e.g. 
Mohr, Fabricant fc Geller 1993 ). Allowing for an ellip- 
ticity introduces two additional free parameters, and the 
information content of some of the images simply is not 
sufficient to provide meaningful constraints. Our goal here 
is to determine whether there is any evidence of regularity 
in intermediate redshift clusters; any mismatch between 
the data and the model (like that in A2390) will serve 
simply as an additional source of scatter in the final rela- 
tion. Therefore, we conservatively present a uniform and 
objective analysis which is appropriate for the majority of 
our sample, rather than varying our analysis from cluster 
to cluster. (An elliptical P model fit to A2390 yields an 
isophotal size which is 8% smaller than that listed in Ta- 
ble 0. This correction makes A2390 more consistent with 
the best fit ST relation, but the correction is unimportant 
compared to the 17% RMS scatter of the sample about 
that best fit relation.) Of course the best approach is 
a non-parametric analysis of the surface brightness maps 
such as that applied to the local sample, but this is not 
possible with the current data. 

We use PROS to convert the detector count rate into 



galactic absorption corrected flux units of ergs s^^cm^^ 
within the cluster rest frame 0.5:2.0 keV band. As for 
the low redshift clusters, we assume a Raymond-Smith 
emission spectrum with the meas ured mean temperature 
Tx and cosmic abundances (Mushotzky & Loewen- 



stein 1997) at the appropriate redshift, and then we cor- 
rect for the (1 + z)* cosmological dimming. The rest 
frame 0.5:2.0 keV central surface brightnesses in units of 
ergs s~^cm~^arcmin~-'^ are also listed for each cluster. 

Angular isophotal sizes 9i are then estimated 
using a fiducial surface brightness of 3.0 x 
lO"""^"* ergs s~"'^cm~^arcmin~^. We determine statistical 
uncertainties in the fit parameters by exploring the like- 
lihood within a grid in the parameters /qj and /?; we 
then use that range to estimate uncertainties in the de- 
rived isophotal size 0/. Figure || contains a plot of the 1, 2 
and 3(7 confidence regions in (3 and Oc for MS1455; lines of 
constant isophotal size 9i are overlaid. Note that lines of 
constant Oj are approx imately parallel to the well know n 
p — 6c correlation (e.g. Mohr, Mathiesen & Evrard 1999), 
minimizing the detrimental effects this correlation has on 
our size measurement uncertainties. 

Measured 9i are then converted into physical isophotal 
sizes Ri using the angular diameter distance cLa'- Ri = 
9 id A- We use the general form 



F 



d-A = 



dz' 
E(z') 



Ho{l + z) 



K = - Qm - 



(14) 



where F{x) = sinh{x), x, and sin{x) in an open, flat 
and closed cosmology. Table |^ contains a list of 9j for 
Ix — 3.0 X 10^^"* ergs/s/cm^/arcmin^ with statistical un- 
certainty estimates and the corresponding Ri for the case 
where VLm = 0.3 & 17a = 0.7. 

Figure ^ contains the intermediate redshift ST relation 
with Ict error bars in both Tx and i?/. The conversion 
from 9i to i?/ assumes = 0.3 and f^A — 0.7 in this 
figure. The best fit relation in this case is 



Ri = (0.67 ±0.04) 



X 



6keV 



1.19±0.21 



Mpc, (15) 



and the RMS scatter in size about this relation is 17%. 
Also plotted is the local ST relation for this same isophote 
(solid line). The uncertainties are derived by bootstrap 
resampling the list of 11 sizes (allowing duplication) and 
refitting the relation. 

There is a suggestion that the intermediate redshift ST 
relation is steeper than the local relation, but the differ- 
ence in slopes is less than la significant. In fact, both 
the zeropoints and slopes of the local and intermediate 
ST relations are statistically consistent when we use this 
particular set of cosmological parameters. With a larger 
sample it will be possible to measure the intermediate red- 
shift ST relation slope more accurately; a comparison of 
the two slopes would be an important diagnostic of un- 
usual structural evolution in clusters. 

Figure || contains a plot of the intermediate redshift ST 
relation assuming = 1- The best fit slope in this case 
is 1.16 ± 0.20, and the zeropoint is 0.59 ± 0.034. This 
zeropoint differs from that of the local ST relation (also 
calculated with fijv/ = 1) at 2.5ct. 
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Table 1 

Cluster Fit Parameters and Isophotal Size 
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We have examined the local and intermediate redshift 
ST relations at four other isophotes: 1.5 x 10"^** cgs, 
2.0 X 10"" cgs, 4.0 X lO-i-* cgs, and 5.0 x IQ'^^ cgs. 
The general conclusions reached above using the isophote 
3.0 X 10"^'' cgs apply equally as well at these other 
isophotes: there is a suggestion that the ST relation at 
intermediate redshift is steeper (at less than 1 a signifi- 
cance) than the local relation, and the zeropoints are in 
good agreement when using the cosmological parameters 
VIm = 0.3 and J^a = 0.7. The fainter the isophote the 
larger the enclosed region. Pushing to fainter isophotes is 
dangerous given the quality of the intermediate redshift 
data, and pushing to brighter isophotes leads to complica- 
tions from central emission excesses and the effects of the 
PSF. Thus, we choose the isophote 3.0 x 10"^'' cgs as a 
compromise between these two competing effects. 

3.3. Cosmological Constraints 

Although the X-ray images of these 11 clusters are poor, 
it is nevertheless interesting to use the local ST relation 
and our evolution model to predict distances to the in- 
termediate redshift sample, thereby constraining cosmo- 
logical parameters. Because the clusters are at a range 
of redshifts, both the scatter of the data about the lo- 
cal ST relation and any systematic offsets from the local 
relation contain cosmological information. Operationally, 
for each set of cosmological parameters VLm and f^A -we 
convert the measured Oj Rj and calculate the of 
the sample about the best fit local ST relation, where the 
best fit local relation is calculated using this same set of 
cosmological parameters. Fig |^ contains contours of Ax^ 
corresponding to 1, 2 and 3 a confidence regions (equiva- 
lent to Ax^ = 2.3, 6.2 and 11.8) in the VLm-^k space. 

The current sample, because of observational uncertain- 
ties, does not provide a strong cosmological constraint. 



The la confidence region is a wide trough similar in char- 
act er to those derived from luminositv distances t o SNe 
la ( Schmidt et al. 1998| , [Pcrlmuttcr et al. 1999|) , with 



the important difference that our sample of 11 intermedi- 
ate redshift distances is not as constraining as the larger 
samples of SNe la distances. Notably, there is enough 
power in the current dataset to exclude the = 1 & 
S^A — model (one of three models marked with stars) 
at >95% confidence. If one considers only models where 
Q,M + f^A = 1, then the preferred model has JIm = 0.09. 
Furthermore, < 0.32 with la confidence and fJ^v/ < 1 
with 3tT confidence. 

One can also use these data to constrain the decelera- 
tion parameter = — f^A/2. The preferred value is 
qo = 0.1ltil;e5 (68% confidence). 

Clearly, larger samples of archival data and the impend- 
ing stream of high quality X-ray images from Chandra 
and XMM will provide tighter constraints. One particular 
concern in this study is that the local and intermediate 
redshift samples were observed with two different instru- 
ments: the ROSAT PSPC and the ROSAT HRI. Any rel- 
ative calibration errors would then serve to bias the pre- 
ferred cosmological parameters. However, the steepness of 
cluster X-ray surface brightness profiles mitigates this po- 
tential problem; typically, the surface brightness falls off 
as / cx 9^^, so even a 10% relative error in the calibration 
of the PSPC and HRI would introduce only a 3% error in 
the angular size of the cluster. 

Although there has been no previous use of cluster scal- 
ing relations to constrain cosmological parameters, there 
have b een three attempts to u s e the ICM mass fraction 
fiCM ( pen 1997| , pooray 1998] , |Rines et al 1999| ); if one 
assumes ficM is constant with redshift, then one can use 
its distance dependence to constrain the distance-rcdshift 
relation. The Pen constraint go = 0.89 + 0.29 is marginally 
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inconsistent with our measurement, but the Pen analysis 
is only appropriate in the case that the ICM traces the 
dark matter, a case that is inco nsistent with both theo- 
retic al and observa tional studies ( David, Jones fc Forman| 
1995|, lEvrard 1997|). 



In the Rines et al. analysis, care is taken to (1) match 
the luminosity ranges of the local and intermediate red- 
shift samples, (2) address the expected radial variation of 
fiCM, and (3) test the consistency of the pubhshed ef- 
fective areas of the two satellites involved (Einstein and 
ASCA). Because ficM increases slowly with radius, it is 
important that it be measured within the same portion 
of the virial region in the local and intermediate redshift 
samples; this complicates the use of ficM measurements 
as cosmological constraints. If one analyzes local and in- 
termediate redshift clusters at a fixed metric aperture, the 
region studied will correspond to a larger fraction of the 
virial region as the redshift increases (see Eqn. ^2|). Simi- 
larly, when measuring ficM within a fixed portion of the 
virial region r^ooi it is important that the redshift evolu- 
tion of the rsoo — be included; it appears that in both 
the Rines et al. and Cooray analyses, no evolution was 
included, so the portion of the virial region used in deter- 
mining ficM increases with redshift. In both cases, these 
biases cause the inferred angular diameter distances to be 
underestimated. 

The size of the bias depends on the cosmological model 
being evaluated (Eqn. n3); using the relation ficM oc 



(^Tir/r^oof -^^ dEvrard 1997|, |Pen 1997D, we estimate the 



bias on at z = 0.35 to be -2% for = 0.3 & Oa = 0.7, 
and -10% for Vl^ = 1 and Oa = 0; these biases are im- 
portant relative to the ^^8% difference in dA{z = 0.35) 
between an JIm = 0.3 and an iljv/ = 0.3 & Ha = 0.7 
cosmology. Correcting for this bias will improve the con- 
sistency between our cosmological constraints and those of 
Rines et al. and Cooray. 

4. CONCLUSIONS 

We present the local X-ray Size-Temperature (ST) re- 
lation for an X-ray flux limited s ample of 45 clusters ob - 
served with the ROSAT PSPC ( |Mohr fc Evrard 19971 ). 



We provide an explanation of this scaling relation in terms 
of underlying scaling relations in the cluster dark matter 
properties. The observed ST relation is slightl y ste eper 
than the self-similar expectation presented in §2.3, but 



is consistent when the modest variation of the ICM mass 
fraction with is taken into account ( Mohr, Mathiescij 
fc Evrard 199^1 ). 

We use our theoretical model for the ST relation to ex- 
plore its evolution with redshift. Interestingly, for the typ- 
ical ICM radial distribution observed in nearby and dis- 
tant clusters, the normalization of the ST relation is not 
expected to evolve. The lack of evolution makes the ST 
relation a plausible source of intermediate redshift angular 
diameter distances. Of course, if cluster structure evolu- 
tion is very different fro m the current theoretical expec- 
tation, our model (§3.1) will underestimate ST relation 
evolution. For example, a shift in the mean ICM mass 
fraction with redshift would bias ST relation distances; to 
date there is no compelling evidence that distant clusters 
have different ICM mass fractions in the mean than nearby 
clusters (e.g. Lewis et al. 199£ , Grego et al. 200C). 



We use ROSAT HRI observations of 11 CNOC clusters 
with measured emission weighted mean ICM temperatures 
Tx to make the first measurements of the intermediate red- 
shift ST relation. Because of the poor image quality, we 
measure the angular isophotal size 9j using the best fit cir- 
cular (3 model, rather than measuringit nonparametrically 
as for the nearby clusters (see Eqn |l|) . By assuming the 
cosmological parameters ^Im = 0.3 « Ha = 0.7, we show 
that the slope and zeropoint of this intermediate redshift 
ST relation is statistically consistent with that of the lo- 
cal ST relation (see Fig In addition, we examine the 
relation for = 1 & ^^A = (see Fig ^), showing that 
although the slope is consistent, the zeropoint is different 
at greater than 3(t significance. 



a 



ST Relation Distances 

1^-3x10"'^ ergs s"' cm"^ arcmin"^ 
Contours correspond to la, 2tJ and 3a 
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Fig. 5. — Cosmological constraints from X-ray ST relation dis- 
tance measurements to 11 intermediate redshift clusters. Contours 
correspond to 1, 2 and 3 a confidence regions (Ax^ = 2.3, 6.2 and 
11.8). The solid line marks spatially fiat models, and the stars mark 
models Qm = 0.3 & Qa = 0.7, Qm = 0.3 and Qm = 1- 

Finally, we use this cluster sample and our ST relation 
evolution model to place cosmological constraints. Given 
the quality of the cluster images and temperature measure- 
ments, it is not surprising that a wide range of cosmolog- 
ical models is consistent with the data. Nevertheless, this 
sample of 11 intermediate redshift distances is sufficient 
to rule out Qj^i ~ 1 with between 2a and 3cr confidence. 
Taken together with ICM mass fraction constraints on the 
cosmologica l matter density parameter Hm < 0.44 at 95% 
confidence (Mohr, Mathiesen & Evrard 199£), the cluster 
ST relation exhibits a slight preference for universes with 
r^A > 0; models with 17 a = are inconsistent with the 
ST relation at between 1 and 2a. When considering only 
models where Hm + Ha = 1, we can rule out f^M = 1 with 
3a confidence. 

With the higher quality X-ray images and ICM temper- 
ature measurements available from Chandra and XMM, 
a significant tightening of these constraints and further 
tests of the underlying evolution model will be possible. 
Comparison of local and distant Mjcm — Tx relations, 
which are more sensitive to cluster evolution, would pro- 
vide important constraints on these models. In addition, 
observations with a new generation of Sunya e v-Zel'dovich 
effect instruments (Carlstrom et al. 199E, Mohr et al 
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1999, Holder et al. 200C ) will allow us to more accu- 
rately constrain the evolution of cluster structure. With 
these future observations of intermediate and high red- 
shift clusters, we plan to continue using the ST relation as 
a tool to provide cosmological constraints independent of 
those derived from re c ent high rcdshift SNe la observations 
( [Schmidt ct al. 1998| , [Perlmutter et al. 1999| ). 



JJM is supported by Chandra Fellowship grant PF8- 
1003, awarded through the Chandra Science Center. The 



Chandra Science Center is operated by the Smithso- 
nian Astrophysical Observatory for NASA under contract 
NAS8-39073. EDR is supported by NASA GSRP Fel- 
lowship NGT5-50173. EE acknowledges support provided 
by the National Science Foundation grant AST 9617145. 
AEE acknowledges support from NSF AST-9803199 and 
NASA NAG5-8458. This research has made use of data 
obtained through the High Energy Astrophysics Science 
Archive Research Center Online Service, provided by the 
NASA/Goddard Space Fhght Center. 



REFERENCES 



Arnaud, M. & Evrard, A.E. 1999, MNRAS, 305, 631 
Bryan, G.L. & Norman, M.L. 1998, ApJ, 496, 80 
Buote, D.A. & Tsai, J.C., 1996, ApJ, 458, 27 

Carlstrom, J.E., Joy, M.K., Grego, L., Holder, G.P., Holzapfel 
W.L., Mohr, J. J., Patel, S. & Reese, E.D. 1999, contributinri tr 
Nobel Symposium on Particle Physics and the Universe, astro- 
ph/c 905255 

Cavaliere, A. & Fusco-Femiano, R. 1978, A&A, 70, 677 

Cavaliere, A., Menci, N. & Tozzi, P. 1998, ApJ, 501, 493 

Cooray, A.R. 1998, A&A, 333, L71 

na.nri T, P Blyz B r. W VrfiloV Bn.tr Arna.iH K A 



1993 , ApJ, 412, 479 

David, L. P., Jones, C & Forman, W. 1995, ApJ, 445, 578 

Edge, A.C., Stewart, G.C., Fabian, A.C. & Arnaud, K.A. 1990, 

MNRAS, 245, 559 
Evrard, A.E., Metzler, C.A. & Navarro, J.F. 1996, ApJ, 469, 494 
Evrard, A.E. 1997, MNRAS, 292, 289 

Fabian, A.C., Crawford, C.S., Edge, A.C. & Mushotzky, R.F. 1994, 

MNRAS, 267, 779 
Fabricant, D., Lecar, M. & Gorenstein, P. 1980, ApJ, 241, 552 
Grego, L., Carlstrom, J.E., Reese, E.D., Holder, G.P., Holzapfel, 

W.L., Joy, M.K., Mohr, J. J., & Patel, S. 2000, ApJ, submitted 
Holder, G.P., Mohr, J. J., Cplstrom, J.E,, Leitch, E.M. & Evrard, 

A.E. 2000, ApJ, in press ( istro-ph/9912364 ) 
Horner, D.J., Mushotzky, R.F. & Schari, C.A 1999, ApJ, 520, 78 
Jones, C. & Forman, W.F. 1984, ApJ, 276, 38, 
Jing, Y.P. & Suto, Y. 1999, ( istro-ph/9909478 ) 
j0rgensen, I., Franx, M. & K aergaard, P. 199b, MNRAS, 280, 167 



Lewis, A.D., EUingson, E., Morris, S.L. & Carlberg, R.G. 1999, ApJ. 
517, 587 

Markevitch, M. 1998, ApJ, 504, 27 
Metzler, C.A. & Evrard, A.E. 1994, ApJ, 437, 564 
Mohr, J. J., Fabricant, D.G. & Geller, M.J. 1993, ApJ, 413, 492 
Mohr, J.J., Evrard, A.E., Fabricant, D.G & Geller, M.J. 1995, ApJ. 
447, 8 

Mohr, J.J. & Evrard, A.E. 1997, ApJ, 491, 13 
Mohr, J. J., Mathiesen, B. & Evrard, A.E. 1999, ApJ, 517, 627 
Mohr, J. J., Carlstrom, J.E., Holder, G.P., Holzapfel, WT,^ w 



_EJVL^k_ai 



ph/9905256| ) 
Moore, B., Q 



E.D. 1999, "VLT Opening Symposium", (astro- 



Lake, G. 1999, 



MNRAS, 310, 1147 (|astro-ph/9903i64) 
Mushotzky, R.F. & Loowenstcm, M. 199 ApJ, 481, L63 
Navarro, J.F., Frenk, C.S. & White, S.D.M. 1997, ApJ, 490, 493 
Pen, U. 1997, New Astronomy, 2, 309 
Perlmutter, S. et al. 1999, ApJ, 517, 565 

Ponman, T.J., Cannon, D.B. & Navarro, J.F. 1999, Nature, 397, 135 

Raymond, J.C. & Smith, B.W. 1977, ApJS, 35, 419 

Reese, E.D., Mohr, J. J., Carlstrom, J.E., Joy, M., Grego, L., Holder, 

G.P., Holzapfe] , W.L., Hughoe, J .P., Patel, S. & Donahue, M. 2000, 

ApJ, 533, 38 ( 3.stro-ph/9912071 ) 
Rines, K., Forman, W., Pen, U., Jones, C, Burg, R. 1999, ApJ, 517, 

70 

Schmidt, B.P. et al. 1998, ApJ, 507, 46 

Yee, H.K.C., EUingson, E., & Carlberg, R.G. 1996, ApJS, 102, 269 



